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We present a continuum theory of self-propelled particles, without alignment interactions, in a 
momentum-conserving solvent. To address phase separation we introduce a scalar concentration 
field (f> with advective-diffusive dynamics. Activity creates a contribution T,ij = —C{(di<t>){dj<j>) — 
(V</>) 2 <5ij / d ) to the deviatoric stress, where is odd under time reversal and d is the number of spatial 
dimensions; this causes an effective interfacial tension contribution that is negative for contractile 
swimmers. We predict that domain growth then ceases at a length scale where diffusive coarsening 
is balanced by active stretching of interfaces, and confirm this numerically. Thus the interplay of 
activity and hydrodynamics is highly nontrivial, even without alignment interactions. 


PACS numbers: 63.50.-x, 63.50.Lm, 45.70.-n, 47.57.E- 

‘Active matter’ means the collective dynamics of self- 
propelled particles at high density. By converting energy 
into motion, such particles violate time-reversal symme¬ 
try (TRS) at the micro-scale. This violation changes the 
structure of coarse-grained equations of motion, allowing 
far-from-equilibrium physics to dominate at large scales 
[Tj. Active matter can be ‘wet’, i.e., coupled in bulk to a 
momentum-conserving solvent, or ‘dry’, i.e., for instance 
in contact with a momentum-absorbing wallj^] ‘Wet’ ac¬ 
tive systems include not only bacterial swarms in a fluid, 
the cytoskeleton of living cells, and biomimetic cell ex¬ 
tracts CHS, but also synthetic colloidal swimmers in a 
fully bulk geometry. Such swimmers may in future of¬ 
fer a toolbox for directed assembly of new materials in 
three dimensions. Many of these artificial colloidal swim¬ 
mers are spherical objects, with asymmetric coatings that 
cause them to move through a bath of fuel and/or under 
illumination with light |5|. 

Of particular importance is the ‘active liquid crystal’ 
(ALC) theory [T] [5]. which starts from a passive fluid 
of rod-like objects [Tj with either a polar order param¬ 
eter P j8], or a nematic one Q M- An active stress is 
then added; this is -CP OP or — CQ, with C odd un¬ 
der time reversal and the dyadic product ®, representing 
the leading-order TRS violation in an orientationally or¬ 
dered medium. This causes new physics such as giant 
number fluctuations nni, and spontaneous flow above an 
activity threshold that vanishes for large systems mm ■ 
This instability depends on whether particles are exten¬ 
sile (pulling fluid inwards equatorially and emitting it 
axially) or contractile (vice versa). Numerical solutions 
HMS1 show spontaneous flows resembling experiments 
on bacterial swarms [2| and on microtubule-based cell 
extracts |3]. 

There is one important effect of activity that ALC 
models do not capture (unless added by hand |il4]): 


1 This wet/dry terminology has become conventional, although 

many ‘dry’ systems are immersed in a fluid. 


motility-induced phase separation (MIPS) [T5, TBj. If 
their propulsion speed falls fast enough with density (e.g., 
due to crowding interactions), even purely repulsive ac¬ 
tive particles phase-separate into dense and dilute do¬ 
mains. MIPS is by now well-established, at least for ‘dry’ 
models such as Brownian dynamics simulations of self- 
propelled colloidal spheres PZH2D]. Unlike the physics 
of ALCs, MIPS does not depend on alignment interac¬ 
tions El. and can be captured at continuum level by a 
scalar concentration field <fi alone, without direct appeal 
to order parameters P or Q [2T1 . This is pertinent 
to spherical colloidal swimmers |5j, which (although not 
entirely devoid of alignment interactions [22]) are not ori¬ 
entationally ordered in the passive limit, contrary to the 
ALC premise. The theory of “scalar active matter", with 
MIPS as its main feature, has been developed so far only 
for dry systems pnaou i23i mi- Following the same path 
as for passive systems manT] culminates in a dynamical 
scalar </> 4 field theory, called ‘Active Model B’ [21] . 

In this Letter we extend this simplest scalar active 
matter model to the wet case. For passive systems it is 
well known how to couple a diffusive, conserved, phase- 
separating order parameter (see passive Model B [251) 
to an isothermal fluid flow; the result is called ‘Model 
H’ [25]. We follow a similar path, but find that activ¬ 
ity deeply alters the relation between the stress and the 
scalar order parameter. At first sight, our final equations 
resemble closely those of passive Model H, in which the 
domain size L(t) scales linearly in time, L ex ert/p , as 
found from a force-balance between interfacial tension a 
and dynamic viscosity 77. On closer inspection though, 
the active equations involve two separate tensions, one 
in the diffusive and one in the mechanical sector. The 
second of these can be (and for pure MIPS actually is) 
negative for contractile particles. As a result, while ex¬ 
tensile systems show relatively normal domain growth, 
contractile ones should cease to coarsen at a certain scale, 
set by a balance between loss of interfacial area through 
diffusive coarsening, and its creation under the action of 
the contractile stress. The latter effect arises because 
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swimming particles are more likely to point normal to 
a fluid-fluid interface than tangential to it (see Fig. 0. 
Within a MIPS context, we also show that this orienta¬ 
tional bias has a purely kinematic explanation: it does 
not require interparticle torques. 
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Figure 1. Schematic diagram showing the flow caused by 
swimmers whose polarization is normal to an interface be¬ 
tween phases. For contractile swimmers this is mechanically 
equivalent to a negative interfacial tension. 


Active Model B: To establish notation and basic con¬ 
cepts, we briefly recapitulate the formulation of [2T1. For 
‘dry’ scalar active matter, the diffusive dynamics of the 
compositional order-parameter field <f>( r, t ) obey0 


0=-V-J, 

J = —V [i T A. , 

fj, = a(j) + be/) 3 4 — kV 2 /> + A(V</>) 2 . 


(1) 

( 2 ) 

(3) 


Here, J is the diffusive current, whose mobility M in 
Eq. § is set to unity, A is a zero-mean, unit-variance 
Gaussian white noise, and a, b, n, and A are constants. 
(Negative a is chosen to ensure phase separation; b and k 
are positive for stability.) This equation set differs from 
the traditional diffusive dynamics of a conserved scalar 
field (see Model B |25| ) solely by addition of a leading- 
order TRS violation, namely the A term. Without this 
term, /x, which resembles a chemical potential, can be 
written as 8T/8cf> with the functional (setting k B T = 1) 

‘ F = /(^ 2 + ^ 4 + ^ ( ' V ^ 2 ) ddr ’ ^ 


where d is the number of spatial dimensions. For ac¬ 
tive systems T is not a genuine free-energy functional, 
and in the simplest MIPS theory, it stems solely from 
the density-dependence of the propulsion speed m- 
Nonetheless, if A happens to vanish, Eqs. 0 ©> coincide 
with Model B, which describes a passive system with free 
energy T [ 25] . In contrast, for nonzero A no functional J- 


In |211 it was shown that Eqs. |IJ-|3]| capture all terms in J 
to order V' ! and o 2 , with one exception, f.r/jV vt, which can be 
removed by introducing a density-dependent mobility; see Ap¬ 
pendix |A] 


exists for which p, = SB/Scfi [21], giving Active Model B. 
Thus only the A term violates TRS in the macroscopic 
equations, even if the physical origins of a, b, and k also 
do so microscopically. 

Active Model H: We now wish to couple the diffusive 
dynamics of (f>{ r, t) to a momentum-conserving solvent 
with fluid velocity v(r ,t). Diffusive dynamics now takes 
place in the frame of the moving fluid so that Eq. 0 
acquires an advective time derivative, 

(j) + v-V/> = —V• J , (5) 

with no change to Eqs. ([2]) and 0. The fluid is incom¬ 
pressible, 


V-v = 0 , (6) 

and of unit mass density. The Navier-Stokes equation for 
momentum conservation then reads 

v + v-Vv = r/V 2 v-Vp + V-E . (7) 

Here, the pressure field p(r, t ) subsumes all isotropic 
stress contributions and enforces Eq. 0. The deviatoric 
stress E is traceless and (without orienting interactions) 
symmetric. 

In passive systems, E can be derived from the free- 
energy functional J- by standard procedures [251 : restor¬ 
ing isotropic terms one finds V-E = —<^V/x which is the 
thermodynamic force density arising from concentration 
gradients m ■ The deviatoric stress can then be written 

s« = -c ((ditm® - (8) 

with ( = k. But if we relax that equality, Eq. 0 remains 
the only deviatoric stress that can be created from </>(r, t) 
to order 0((/ 2 , V 3 )0 and is hence the sole leading-order 
deviatoric stress contribution for scalar active matterQ 
This neither contradicts, nor depends on, recent analyses 
of the pressure and/or stress in specific active models 123- 

m- 

Our Active Model H comprises Eqs. 0, 0, and 0- 
0. Regardless of the underlying cause of phase separa¬ 
tion (conventional attractions and/or MIPS in any com¬ 
bination), only two activity-dependent parameters dis¬ 
tinguish the dynamics of Active Model H from a passive 
Model H with the free-energy functional J- given in Eq. 
0 > namely A and 


3 This can be established by constructing all traceless and sym¬ 
metric second-rank tensors T to that order, and showing that 
each of their divergences varies as VT up to the gradient of a 
scalar field (i.e., a pressure gradient); see Appendix |B| 

4 The parameter £ = Cr + Cd contains now a reversible contri¬ 
bution (even under time reversal) and a dissipative one 
(odd under time reversal) that stems solely from activity. If lin¬ 
ear irreversible thermodynamics holds, we have = k. This is, 
however, not the case for systems that are sufficiently far from 
equilibrium m ■ 
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Our previous work on Active Model B shows A to have 
benign effects on dynamics: its main effect is to shift the 
coexistence condition between phases which now obey 
an ‘uncommon tangent’ construction PH]. Once this 
shift is allowed for, standard diffusive dynamics (e.g., 
Ostwald ripening mi) holds qualitatively; although sim¬ 
ulations show deviations from the expected power law 
L oc t 1//3 , these likely represent a slow crossover to that 
asymptote m- Below we thus focus mainly on the role 
of £. 



Figure 2. Results for the domain size L{t) in extensile (£ > 0) 
and contractile (£ < 0) two-dimensional scalar active fluids. 


Physics of the stress term: This can be addressed in 
two contrasting limits. Suppose first we have a sys¬ 
tem with strong interparticle attractions, which would 
phase separate without activity. Then the thermody¬ 
namic, zeroth-order result is £ = k as mentioned previ¬ 
ously. Perturbative activity will change this only slightly. 
TRS is broken, but as with A the effects on phase order¬ 
ing dynamics should be benign. Specifically, we can use 
standard power-counting arguments E7IEU to estimate 
the dynamics of the domain size L(t). Setting A = 0 for 
simplicity and ignoring inertia (which leads to a third 
regime at very late times j3H 03), we find that (see Ap¬ 
pendix [C]) 

L ~ aa/L 2 + f3a/p . (9) 

Here, a and (3 are dimensionless constants; <7 = 
(— 8na 3 /9b 2 ) 1 / 2 is the interfacial tension in </> 4 theory. 
When activity is nonzero but small, C, ^ n and the ten¬ 
sion a differs slightly from er. The two contributions to L 
in Eq. ^ are related to gradients in p, and thus Laplace 
pressure n. The first is diffusive (so L is a flux propor¬ 
tional to V^t ~ Il/L ~ er /L 2 , where er is unperturbed 
since A = 0); £ does not enter here. The second is hy¬ 
drodynamic (balancing a viscous stress r/V 2 v ~ r\L/L 2 
against V • S ~ d/L 2 ). When a = er, Eq. Q cap¬ 
tures the well-studied crossover from Model B behavior, 
L oc (crt) 1 / 3 , to viscous hydrodynamic (VH) coarsening, 
with L cx (tt/r), at a length scale L x oc r\ x ! 2 . Despite loss 


of TRS, a perturbative shift in a is unlikely to change 
this outcome, since Eq. (J9j) rests on little more than di¬ 
mensional analysis EH- 

The opposite limit is that of pure MIPS, where there 
are no attractions between particles and every coefficient 
in Active Model H is set by far-from-equilibrium physics 
(primarily, the dependence of a particle’s swim speed on 
the local density mi isi). In a TRS system the form of 
J- stems from a Hamiltonian that determines the diffu¬ 
sive (from V/r) and mechanical (from £) currents in a 
thermodynamically consistent way. However, as already 
mentioned, in MIPS, T is merely a mathematical con¬ 
struct; there is no Hamiltonian (even for A = 0) and 
hence no link between n and £. Indeed, whereas k is 
always positive, £ can have either sign as we now show. 

The argument is based on simple kinematics as devel¬ 
oped in m ■ For a system of swimmers with a propulsion 
speed w(p) that depends on position (or particle density 
p) but not on swimming direction u, the first and sec¬ 
ond orientational moments of the distribution function 
^(r, u) obey |Hj 


P = —rV(wp) , 

Q = vm> )- 


( 10 ) 

(ii) 


Here, r is the orientational relaxation time; the par¬ 
ticle density p is the zeroth orientational moment of 
\F(r, u); and ST denotes the symmetric traceless compo¬ 
nent. Note that p and 4> are related linearly; specifically if 
one sets a = b = 1, one has £> = (2p — pH — Pl)/{ph — Pl), 
where pp and p p denote the low and high densities of co¬ 
existing phases, respectively m■ Equation is purely 
kinematic in origin: wherever wp has a gradient, there 
are more particles pointing (i.e., moving) down this gra¬ 
dient than up it, causing nonzero P [24] . 

For self-propelled particles, the leading-order mechan¬ 
ical stress is, as for ALCs, caused by their exerting force 
dipoles on the fluid [I]. In general we can write this as 
£ = —£pP®P/p—£ qQ where £p and £q are activity pa¬ 
rameters. However, from Eq. (Ill it follows that if w is a 


function of density only (as assumed to this order in the¬ 
ories of MIPS [T9] ), V-Q is a pure pressure gradient and 
thus ignorable. Expanding Eq. (101 as P = —T(wp)'S7p 
with prime denoting a p derivative, we recover Eq. (|8j, 
where 


C(0) 


Cp / r(wp)' \ 

p \ V ) 


( 12 ) 


The sign of £ then depends only on whether swimmers are 
extensile (£p > 0) or contractile (£p < 0). Time-reversal 
interchanges these two cases, so the active contribution 
to £ in Eq. ([8]) is odd under it, unlike any passive part(^] 


5 Note that Eq. § satisfies the general criterion for an admissible 
active stress presented in m- 
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Finally, although £p vanishes for spherical swimmers at 
low density E3 it does not do so in genera© 

In keeping with our earlier discussion of Eq. ([ 8 ]) we 
now suppress the (/-dependence of £. However, one could 
alternatively retain £(£>) in Eq. as part of a ‘best- 
fit’ procedure to a more detailed kinetic theory of MIPS 
(see Appendix]©; we have checked numerically that the 
results are broadly similar to those with constant £ re¬ 
ported in what follows. 

Next, we assume that local diffusive relaxation nor¬ 
mal to interfaces ensures that their local structure is only 
weakly perturbed by fluid motion, just as applies in pas¬ 
sive Model H E2- We thereby nonperturbatively recover 
Eq-©, now with a = (a/n (see Appendix |C]). So long as 
swimmers are extensile, the two interfacial tensions in Eq. 
© are positive, and the diffusive ( L oc (ert) 1 / 3 ) and VH 
f L oc crt/rf) regimes both remain intact. The prefactor of 
the VH coarsening is however shifted (possibly by a large 
amount), as is the crossover length L x oc (rja/a) 1 ^ 2 . 

In contrast, very different physics now arises for con¬ 
tractile swimmers. Here, the interfacial tension a in the 
diffusive sector is positive as usual, but the mechanical 
one, a, is negative. Equation © still makes sense, but 
instead of a crossover from diffusive to VH coarsening 
there is a balance point Lb oc ( 770 -/|<t| ) x / 2 where diffusive 
coarsening is negated by VH anti-coarsening. Thus we 
predict that growth of domains will cease at this scale, to 
be replaced by a dynamic equilibrium where the diffusive 
shrinkage of interfacial area is balanced by its production 
under the action of the contractile stress. This stretch¬ 
ing of interfaces can be understood as a combination of 
the kinematics - causing swimmers at an interface to be 
aligned normal to it on average - with the flow pattern 
around a contractile swimmer, which pulls fluid inwards 
axially and expels it around the equator. For a swimmer 
normal to the interface (with either polarity) the latter 
is equivalent to a negative interfacial tension (see Fig. 
©. This argument also shows why any (/-dependence of 
£ (unless it changes sign) should not alter things quali¬ 
tatively. Negative active tension has also been suggested 
in [33] by an argument that appears unrelated, in that it 
applies equally to contractile and extensile cases. 

Numerical results: To test these ideas we have solved 
Active Model H numerically on a square lattice of size 
l x = l y = 256 by using a noise-free (i.e., A = 0) hybrid 
Lattice Boltzmann (LB) scheme [3S]. We initialized the 
system in a mixed state with a small initial uniformly 
distributed noise with —0.1 < <j>( r, 0) < 0.1. In our sim¬ 
ulations we set a = —0.004, b = 0.004, k = 0.006, A = 0, 


For an isolated isometric swimmer with orientation u the only 
tensor from which to construct a stress contribution is u(8>u, 
whose mesoscopic average at low but finite density gives Q but 
not P(x)P w- At high density, where swimmers interact strongly 
and correlations develop between orientations and interparticle 
coordinates, the reasons given in m for the vanishing of the 
P 6Cj P term no longer apply. 



Figure 3. Time series showing growth and arrest of domain 
structure in a two-dimensional contractile scalar active fluid 
with £ = —0.001 and box size 256 x 256. 


and 77 = 1.67, and varied £. We checked that setting 
A 7 ^ 0 does not qualitatively affect our conclusions. 

These simulations confirm our expectation that con¬ 
tractile (but not extensile) activity should lead to the 
arrest of MIPS (see supplemental movies). This can be 
seen in Fig. 2] which shows Lit) for active systems with 
various £. Here, L(t) is defined as usual through mo¬ 
ments of the structure factor [3?|. Snapshots of the dy¬ 
namics corresponding to an arrested state are shown in 
Fig-! and in supplemental movies; although the length 
scale is now fixed, the structure is highly dynamic. As 
£ increases towards zero, the saturation length grows 
and for —10~ 4 < £ < 0 coarsening does not arrest 
within the time window of our simulations. We spec¬ 
ulate that it would do so eventually (at a length scale 
diverging as £ -7 0 _ ) although our data do not rule out 
a finite negative threshold above which arrest never oc¬ 
curs. The saturation length Lb was argued above to 
obey Lb oc (r]a/\a\) 1 ^ 2 . At fixed 77 and k this gives 
Lb oc |CI~ 1/2 - O ur data for Lb{() instead approach a 
weaker power law at low activity and a plateau at larger 
values (see Fig. [4] in Appendix |E| . The cause of this dis¬ 
crepancy is unclear, but might suggest that a force bal¬ 
ance different from Eq. © prevails in the highly dynamic 
states observed in strongly contractile systems. 

The case £ = 0 equates to passive Model B (no coupling 
to fluid motion), whereas passive Model H is attained 
for £ = k. The L{t) curve is as expected in each case 
with final power laws close to 1/3 and 1, respectively. 
For passive Model H the crossover in L(t) from early- 
time diffusion to late-time fluid flow is sigmoidal, and 
the negatively curved part has often been interpreted as 
a sublinear power law |HJ. Whatever the precise inter¬ 
pretation is, extensile systems with 0 < £ < k interpolate 
smoothly between the two pseudo-passive limiting cases. 
Sigmoidal curves continue to be seen for systems with 
£ > k, in which the tension in the fluid sector is higher 
than the one driving diffusion. The data does not exclude 
an eventual power that exceeds unity, but more likely re¬ 
flects a prolonged sigmoidal crossover between diffusion 
and a linear growth with a large slope L oc cr/ 77 . 
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In summary, we have constructed Active Model H, 
a minimal model for scalar active matter coupled to a 
momentum-conserving solvent. TRS is broken by two ef¬ 
fects. One, encoded by A, causes shifts in the densities 
at which phases can coexist in diffusive equilibrium EH- 
The second effect is new, and amounts to a mismatch 
between the interfacial tension a that drives fluid mo¬ 
tion and the tension a that drives diffusive fluxes. For 
extensile particles (or weakly contractile ones in cases 
where phase separation is driven by attractive interac¬ 
tions rather than being purely motility-induced) both 
tensions are positive, and while their inequality violates 
TRS, its effects on coarsening dynamics are expected to 
be quantitative, not qualitative. The opposite is true for 
contractile particles which become aligned normal to the 
interface between phases, creating a flow pattern that 
stretches it (cr < 0). Balanced by diffusion, this effect 
can cause the domain size L(t) to saturate. 

Finally, it is tempting to associate our prediction of 
arrest at a finite length scale Lb in wet contractile scalar 


active fluids with observations of finite cluster forma¬ 
tion, rather than full phase separation, in synthetic col¬ 
loidal swimmers undergoing MIPS [5]. However, there 
are two objections to this. First, we do not know which, 
if any, of these systems are contractile. Second, most ob¬ 
servations of cluster phases are in ‘nearly dry’ systems: 
clusters are found within two-dimensional layers close to 
a momentum-absorbing boundary [5]. Nonetheless, our 
prediction is that arrested separation should be generic 
in contractile wet systems undergoing MIPS. We look 
forward to future experimental tests of this prediction. 
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Appendix A: Form of the diffusive current 


Appendix B: Form of the stress tensor 

To second order in </> and third order in V the diver¬ 
gence of a traceless and symmetric second-rank tensor 
Tij can always be written as 

djTn = codjdd^id^)) + (V0) 2 , V 2 </>) (Bl) 

with a constant cq and a scalar function (V0) 2 , V 2 </>) 
that depends on <j>, (S7(f>) 2 , and V 2 </>. This means that if 
we neglect multiples of the identity matrix, the stress 
tensor E^ can be written as 


To second order in /> and third order in V the general 
form of the diffusive current J in Eq. ([2| is 

J = -V W</») 2 ) +e</>V 3 <)>+A+0(</> 3 , V 4 ) (Al) 

with a functional 

J7 = J(fM + ^r (' v ^ 2 ) d " r > ( A2 ) 

where /(</>) and k(4>) = kq + <j> are scalar functions, 

Ao, e, kq, and k\ are constants, and A is - in accordance 
with the fluctuation-dissipation theorem - a zero-mean, 
unit-variance Gaussian white noise. Equation (All is 
equivalent to 

J = -v(/'fo)-«($vV + A(V0) 2 ) 

+ e^+ A + 0(</> 3 ,V 4 ) 


with A = Ao — «i/2. We can rewrite Eq. (A31 as 

J = -M($(v(/'(</>) - £(<£)V 2 0 + \mv</>)*) 

- A'O/OV^V^) 2 ) + A + 0(</> 3 , V 4 ) 

with the density-dependent mobility 

M(<j>) = («(«£) + e<j)) e ^ Kl+t 

and the scalar functions 

'f'W 


f'W = 
k{ 4>) = 
H<t>) = 


J M{<t>) 

k( 4>) + e 

M{4) 

A 


# , 


M{<j>) ' 


(A4) 

(A5) 

(A 6 ) 

(A7) 

(AS) 


When we approximate M{(f >) in Eq. (A4) by a constant 
and discard all contributions that are of third or higher 
order in <fi , we obtain 

J = -V(/'(0) - + A(V^) 2 ) + A (A9) 

with k{<f>) = k 0 + («! + e)cj>. Choosing /(</>) = acj) 2 / 2 + 
6</> 4 /4 and approximating the function k(4>) by a con¬ 
stant k finally leads to Active Model B [see Eqs. @-©l- 
Notice that A is now in accordance with the fluctuation- 
dissipation theorem although the same does not apply to 
Eq. (|A4|. 


= c 0 (dW(dW +OW, V 3 ) . (B2) 


Appendix C: Dynamics of the domain size 


Here we present the derivation of Eq. ([9]), which de¬ 
scribes the dynamics of the domain size L(t ) ~ by standard 
power-counting arguments. For this purpose, we consider 
the domain size in the diffusive and in the hydrodynamic 
regime, but neglect the inertial regime. 

In the diffusive regime, the dynamics of the domain 
size follows from the steady-state version of Eq. © to¬ 
gether with Eq. where the noise A can be neglected: 
v-V</> = V 2 /r. To estimate the dynamics of L(t), we 
replace v —> L, V —> 1/L, and <^ —» 1. We further 
set A = 0, which allows to replace /i -A tr/L, where 
a = {-8na 3 /9b 2 ) 1 / 2 is the interfacial tension in Model 
B EH- These replacements lead to the dynamics in the 
diffusive regime: L oc a/L 2 . 

In the hydrodynamic regime, in contrast, the dynamics 
of the domain size follows from the inertia-free version of 
Eq. 0 : ? 7 V 2 v = 'S/p — V-E. We now replace v —>• L, 
V -A 1/L, and pd l3 — E, (? -A a/L [33], where the tension 
a is different from a due to activity. This leads to the 
dynamics in the hydrodynamic regime: L oc a / 77 . 

The full dynamics of the domain size L{t) is estimated 
by a linear interpolation between the diffusive and the 
hydrodynamic regimes, 

L ~ aa/L 2 + fid/r) , (Cl) 


where a and (3 are dimensionless constants. If we assume 
that local diffusive relaxation normal to interfaces en¬ 
sures that their local structure is only weakly perturbed 
by fluid motion, as is the case in passive Model H [27] . 
a in the first term on the right-hand side of Eq. (Cl) is 
not affected by cr in the second term deviating from er. 
Since £ from Eq. ([ 8 | equals k and ef in Eq. (Cl I equals a 
in passive Model H, we can now approximate a = (/ 0 /n. 


Appendix D: Explicit approximation of C(</>) 

An estimate of £(</>) is found by assuming a pure MIPS 
mechanism, and inferring from the Model H free energy 
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an approximation for the density-dependent swim speed 
w(p) (or equivalently w (<"/>)) that appears in the micro¬ 
scopic expression for (((j>) [see Eq. 01- 

First, we use the linear transform stated further above 
to relate the order parameter <J) in Model H to the particle 
density p in a system of self-propelled particles undergo¬ 
ing MIPS. This transform can be rewritten 


P = 



(Dl) 


with the central density p = (pn + Pl)/2, where p l and 
Ph are the low and high coexisting values of p in the 
phase-separated state, respectively. When A = 0, as con¬ 
sidered here, these are found from a common tangent 
construction on the effective free-energy density [15l [24] 


The constant a can then be expressed as 


- _ PH - PL A dln(w) 

2 P \ dln(p) 



(D4) 


This completes the specification of w(p). The result can 
be used in Eq. ( |T2| if one wants to replace the simplest 
version of Active Model H (which has constant () with 
an estimated (((/>) found by assuming that the quartic ef¬ 
fective free energy arises solely from a MIPS mechanism. 

If preferred, Eq. (D3) can be converted into a in¬ 
dependent propulsion speed u>(</>) using the linear trans¬ 
form (Dl I: 


w{4>) 


2pw{Q)e “ Mo 
2 P + Hph - Pl) ' 


(D5) 


f P = P(ln(p) ~ 1) + / In (w(u)) d u . 


(D2) 


Appendix E: Saturation length 


The density-dependent propulsion speed w(p) of the ac¬ 
tive particles, which is needed to fully specify P and Q, 
can be found in terms of (j> by relating the two nonequi¬ 
librium chemical potentials p p = f' p {p) = ln(p)+ln(ti>(p)) 
and p 0 = = a<j> + be/) 3 , respectively expressed 

in the p and (f) representations. In the following, we 
choose a = —b= —1 for clarity, but the calculations can 
straightforwardly be extended towards general constants 
a and b. Since p and </> are related linearly, so should be 
their chemical potentials: we write p p = —apo + b with 
constants a and 6^] This relation is equivalent to 


g-®A*o 

w(p) =pw{p)—— (D3) 


with pw(p) = exp (b), where the constant w(p) can be 
estimated directly, or absorbed into £p in Eq. 


The constant a should be chosen to match as closely 
as possible the bulk chemical potential contribution po = 
— <j> + (j) 3 in the (f > 4 theory to the non-polynomial form 
found from Eq. (D2) via Eq. (Dl I. This matching cannot 
be perfect and the best-fit criterion is therefore somewhat 
subjective. We choose to match the logarithmic deriva¬ 
tive of the swim speed at the central density p of the 
phase coexistence, dln(w)/dln(p)| p= p = pw'(p)/w(p). 


Figure |4] shows the saturation length Lb{ C) for con¬ 
tractile active particles (i.e., £ < 0). The data from our 
Lattice Boltzmann simulations show a weak power law 
Lb = 8|C| -0 ' 105 at low activity (0.002 < |£| < 0.01) and 
a plateau for larger activity (0.03 < |£|). The exponent 
« —0.1 of the power law is larger than predicted by our 
simple theory, which yields an exponent of —0.5. 


'eT 


-7 -6 -5 -4 -3 

MICI) 

Figure 4. Lattice Boltzmann simulation data (red disks) and 
power-law fit curve (blue straight line) for the saturation 
length Ls as a function of |£| for contractile activity (£ < 0 ). 



When T is the free-energy functional corresponding to the 
chemical potential fi p = <5L'[p]/(5p and p = 70+71 </> and (j) 
are related linearly, the chemical potentials p p = SJ^lpj/Sp = 


(I /71 )ST[4>]/S(/) = (l/ 7 i )//0 and = 6T[<j>]/6(f> are propor¬ 
tional. Taking into account that chemical potentials are defined 
only up to a constant offset then yields a linear relation between 
p, p and P(p- 














